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PROG RES  DANS  LA  CONNAISSANCE  DE  LA  TURBULENCE  ATMOSPHERIQUE 
par  C.  Coupry 

Office  National  d' Etudes  et  de  Recherches  A£rospatiales  (ONERA) 
92322  ChStlllon  Cedex  (France) 


RESUME 


L'expos£  tralte  de  mesures  fines  de  turbulence  effectu£es  en  vol,  sur  avions 
sp£cialement  £quip£s,  alnsl  que  des  statistiques  de  masses  d£dultes  des  accelerations 
mesur£es  sur  avions  commerclaux  au  cours  de  vols  r£guliers.  Dans  ce  dernier  cas,  on 
effectuera  une  analyse  critique  des  diverses  m£thodes  de  restitution  de  la  turbulence 
3  partir  des  acc£l£rations,  et  on  sugg£rera  un  certain  nombre  d'actions. 

11  est  montr£  que  les  paquets  de  turbulence  ob£issent  3  une  distribution  de 
Poisson. 


PROCRESS  IN  THE  ANALYSIS  OF  ATMOSPHERIC  TURBULENCE 


ABSTRACT 

The  report  deals  with  measurement  of  time  histories  of  turbulence,  made  by 
specially  equipped  aircraft,  as  well  as  with  statistics  of  events  deduced  from  a 
large  number  of  commercial  flights.  A  critical  analysis  of  the  methods  used  to  reduce 
the  flight  data  is  developed,  and  proposals  for  Improvements  are  suggested. 

It  is  shown  that  patches  of  turbulence  for  which  a  given  level  is  exceeded  are 
distributed  following  a  Poisson  distribution. 


I  -  INTRODUCTION 

Apr? s  une  p£riode  d' Intense  actlvlt£  dans  les  ann£es  60,  les  recherches  sur  la  turbulence  atmosphi- 
rlque  et  son  influence  sur  le  pilotage  et  les  charges  ont  £t£  ralenties  —slnon  abandonees-  dans  beaucoup 
de  pays  :  les  m£thodes  de  certification  £talent  bien  accept£es,  et  les  nouveaux  avions  ne  semblalent  pas 
poser  de  probldmes  particulars. 

Le  regain  d'lnt£r£t  pour  le  probldme  depuls  le  d£but  des  ann£es  80  tient  3  notre  avis  aux  raisons 
suivantes  : 

-  Projets  d'avlona  de  transport  civile  ou  milltaires  d'allongement  de  plus  en  plus  grand,  et  de  struc¬ 
tures  tr4s  souples, 

-  Generalisation  de  l’utlllsatlon  du  contrSle  actlf  g£n£ralls£,  en  particuller  pour  le  contrSle  des  char¬ 
ges  de  manoeuvre  et  de  turbulence, 

-  Missions  de  p£n£tratlon  d'avlons  milltaires  3  tr?a  basse  altitude  et  3  grande  vltesse,  qul  n£cessltent 
un  contrSle  de  vol  compatible  avec  le  sulvl  de  terrain, 

-  Incoherence,  dans  certains  cas,  des  conditions  de  certification  d£dultes  des  concepts  de  rafale  lsol£e 
ou  de  turbulence  continue,  etc... 

Pendant  la  mSme  p£rlode,  des  progris  considerables  ont  £t£  effectu£s  dans  le  domalne  des  m£thodes  de 
mesure,  dans  les  posslblllt£s  des  ordinateurs  et  dans  le  traltement  du  signal.  L'exlstence  de  besolns,  et 
de  moyens  pour  les  aatlsfalre  au  molns  partiellement?  sont  3  l'orlglne  d'une  actlvlt£  de  nouveau  en 
expansion  depuls  quelques  ann£es. 

Les  prlnclpaux  th?mes  d'actlvlt£  d£velopp£s  actuellement  en  France  portent  : 

-  sur  la  mesure  fine  de  la  turbulence  3  partir  d'avlons  ap£clalement  £qulp£s, 

-  sur  la  crSatlon  d’une  banque  de  donn£es  d£dultes  des  acc£l£rations  mesur£es  au  cours  des  vols  r£gullers 
de  la  flotte  de  British  Airways, 

-  sur  la  d£termlnatlon  expirlmentale  des  efforts  dus  3  la  turbulence,  en  longitudinal  et  lat€ral,  par 
essals  sur  maquettes  dans  une  soufflerie  tranasonlque  £qulp£e  d'un  g£n£rateur  de  rafales, 

-  sur  la  pr£vlslon  par  calcul  des  r£ponses  de  1 'avion,  en  pr£sence  d'une  turbulence  cyllndrlque  ou  iso- 
trope, 

-  sur  1' identification  en  lnstationnalre,  en  vue  du  contrOle  actlf,  de  dlff£rents  types  de  gouvernes 
(ailerons,  flapperons,  spoilers), 

-  sur  le  dlveloppement  de  lots  de  contrSle. 

Nous  lialterons  l'expos£,  dans  cet  ACARD- REPORT,  aux  deux  premiers  points,  c'est-3-dlre  aux  mesures 
fines  et  aux  mesures  de  masse  de  la  turbulence  atmoaphSrique. 
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11  -  MESURES  FINES  DE  LA  TURBULENCE  ATMOSPHERIQUE  I 

I 

A  la  fin  da*  annCes  60,  plualeura  caapagnea  de  aesure  dlrecte  de  la  turbulence  on  (tt  effectives  en  l 

France  aur  Mirage  3A,  Pari*  2,  Mystdre  20  et  Concorde,  2  partlr  de  girouettea,  d'accSICroadtrea  et  de 
gyroaitrea  ;  dana  le  courant  dee  annSes  70,  l'ONERA  a  equip!  lea  prototype*  dea  Airbus  de  systdaes 
aeablables,  qul  ont  8tS  occaalonnelleaent  utilises.  L* experience  aontra  que  lea  aeaurea  effectives  3 
l’dpoque  n’Stalent  aalnes  que  Jusqu'2  une  frequence  d'envlron  5  Hz. 

i 

Lea  aotivatlona  pour  une  reprlae  de  cea  aeaurea  fine*  de  la  turbulence  sont  eaaentlelleaent  :  1 

-  la  recherche  d' inf ornat Iona  aalnes  Jusqu'2  dea  frequences  de  10  Hz,  en  vue  du  contrSle  actlf  en  boucle 
ouverte  d'avlona  allltalres  lora  de  p6n!tratlons  2  basse  altitude, 

-  l'6tude  de  la  coherence  de  la  turbulence  en  envergure  pour  appr!cler  lea  posslbllltCs  du  contrfile  actlf 
dea  charges  en  boucle  ouverte  aur  avion  civil  de  grand  allongement. 

Lea  essala  ont  !t!  effectiVs  au  Centre  d’Esaala  en  Vol  de  Istrea  respectiveaent  aur  Mirage  3B  et 
Nord  260. 

II. 1  -  Principe  dea  aeaurea  : 

La  coaposante  vertlcale  W  de  la  turbulence  eat  explicltee,  dana  l'hypothdse  dea  petite  aouveaents, 
par  la  foraule  clnEaatlque  clasalque  : 

w  *  vs  e  _  va  a  +  z 

o0  01  est  l'lncldence,  mesurde  par  une  glrouette, Vo  la  vltease  aCrodynaalque,  obtenue  par  aesure andaood- 
trique, Vj  la  vltease  vrale  par  rapport  au  aol  et  Z  la  vltease  vertlcale  de  l'avlon. 

En  pratique,  on  aesure  l'accei!ratlon  Z  et  la  vltease  de  tangage  6  ,  et  la  turbulence  eat  restitute 
par  la  foraule  : 

w«  Az+v,  e)dr_v0a 

-4 

Af In  d'Cvlter  les  derives  d' Integration,  le  r!sultat  est  flltrC  par  un  flltre  passe— hsut  (calC  2  la 
frequence  de  0 ,25  Hz) . 

L' Installation  dea  capteura  de  aesure  aur  Mirage  3  eat  prCsentCe  dans  la  figure  1  ;  pour  le  Nord 
260,  lea  aeaurea  de  turbulence  sont  faltea  en  trols  points  (nez  de  l'avlon  et  extr6alt6s  d'alles)  ;  le 
sch!aa  d* Installation  eat  donnC  dans  la  figure  2. 

Lea  progrds  vers  une  aellleure  aesure  de  la  turbulence  ont  CtC  obtonua  : 

-  par  un  Etalonnage  soigrV  de  la  glrouette  en  aoufflerle. 


-  par  une  correction  de  la  fonctlon  de  transfert  du  gyroadtre,  dont  la  bande  passante  !tait  trop  falble, 

-  par  une  prise  en  coapte  prCclae  de  la  aodlflcatlon  de  l'lncldence  locale  due  2  la  presence  du 
fuselage. 

Cette  dernldre  correction  a  !t!  effectuCe  lor*  d'essai*  en  vol  en  air  tris  calae,  pour  lesquels  on 
s'attendalt  2  aeaurer  une  turbulence  nulle  ;  dans  cea  condition*,  le  pilots  aolllcltalt  l'avlon  2  dlf- 
fdrentes  frequences,  et  on  aodlflalt  le  gain  de  la  glrouette  jusqu'2  c e  que  la  mesure  de  turbulence 
s'!vanoulsse  pratlqueaient.  Le  coefficient  de  correction  de  glrouette  (1,79  pour  le  Mirage  3)  a  6t!  d!ter- 
■lnf  dans  cea  conditions. 

Les  figures  3  et  A  reprtaentent  le  rtsultat  de  cet  (talonnage  de  glrouette  en  air  calae,  aur 
Mirage  3,  pour  dea  solllcltat ions  du  pilots  2  deux  frequences  dlffdrentes. 

II. 2  -  Rtsultat a  obtenus  aur  Mirage  3  : 

Plualeura  essala  en  vol  ont  ttf  effectuis  aur  Mirage  3  dans  les  conditions  sulvantes  : 

Z-SOOfr  •  Va  m  500  kmf* 

Z.  500  ft  ■  V0  .  300  kmti 

Au  cours  de,  ces  essats,  on  a  aesurt  slaultantaent  les  paraattres  permettant  la  restitution  de  la 
turbulence  (tt,  •  et  Z  ),  alnsl  que  l'acctlfratlon  au  centre  de  gravltt  et  2  la  polnte  avant.  Lea  rtaul- 
tats  sulvants  ont  ttt  obtenus,  : 

"  Spectres  directs  de  turbulence  :  un  spectre  est  prfaentt  dans  la  figure  5  en  tchelle  log-log,  et 
coaparf  au  spectre  thforlque  de  Karaan  (en  polntlllt)  ;  on  note  un  accord  excellent  entre  thtorie  et 
exptrlence,  qul  seable  assurer  que  les  aesure*  de  turbulence  sont  aalntenant  salnes  jusqu'2  plus  de 
10  Hz.  (Les  rlsultats  sont  faussts  par  la  presence  du  flltre  passe-haut  pour  les  frequences  lnftrleures  2 
0,32  Hz.) 

-  Spectres  directs  d'accClCratlons  (en  particular  acciKratlon  de  la  polnte  avant)  ;  un  de  ces  spectres, 
pr<senti  dans  la  figure  6,  aontre  netteaent,  entre  8  et  12  Hz,  les  rCponses  de  structure  dues  2  la 
flexion  fondaaentale  de  la  vollure  et  2  la  flexion  vertlcale  du  fuselage. 


1-3 


-  Spectres  crolsSs  entre  _ accelerations  et  turbulence 
sortie  (une  acceleration  Z 
par  la  formule  classlque  : 


connalssant  1'entrSe  (la  turbulence  H)  et  la 

Twi 


t  5WZ(t) 

'wZ  \  t>-  -Z - 

5wwll ) 


oO  Swz  ( t )  est  la  densite  spectrale  crois€e  entre  W  et  Z  et  la  densite  spectrale  dlrecte  de  la 
turbulence.  La  figure  7  illustre  ce  rEsultat  et  presente,  en  fonction  de  la  frequence,  les  parties 
rSelles  et  lmaginaires  de  cette  fonction  de  transfert. 


-  Finalement,  on  a  profits  de  ces  mesures  de  rSponse  de  1* avion  4  une  turbulence  connue  pour  vallder  les 
codes  de  calcul  de  cette  rSponse  :  un  essai  de  vibration  au  sol  a  permis  de  determiner  avec  precision  les 
caracteristlques  modales  de  l'avion  essay£,  dans  sa  configuration  de  vols  ;  le  calcul  des  forces 
aerodynamlques  lnstatlonnaires  a  6t£  fait  par  une  mSthode  de  doublets.  La  figure  8  presente  la 
comparalson  "calculs-essais"  qul  est  tr4s  satisf alsante  dans  la  bande  de  frequence  explorSe. 

11.3  -  RSsultats  obtenus  sur  Nord  260  : 


Le  contrSle  des  charges  de  manoeuvres  et  de  turbulence  est  un  impSratlf  auquel  dolvent  falre  face 
d£s  malntenant  les  bureaux  d’etudes  qul  conqolvent  les  nouveaux  avlons  de  transport,  tous  dot6s  de  tr4s 
grands  allongements.  L'objectif  des  essais  en  vol  sur  Nord  260  Stait  d'apprScier  dans  quelle  bande  de 
frequences  un  contrSle  en  boucle  ouverte  Stait  envlsageable,  c'est-4-dire  dans  quelles  conditions  une 
information  locale  sur  la  turbulence  Stait  representative  de  la  rafale  rencontrSe  par  l’avion.  On  a  done 
cherchS  4  mesurer  la  turbulence  en  plusieurs  points  de  l’apparell  (perches  d'extrSmltS  d'aile  et  nez  de 
l'avion)  et  4  determiner  la  coherence  des  mesures. 

Plusieurs  essais  ont  StS  effectuSs,  3  150  kmts  et  180  kmts,  pour  une  altitude  de  500  ff  .  ;  au  cours 
de  ces  essais,  le  mesure  centrale  n'a  pas  fonctlonnS. 

On  prSsente  lcl  les  rSsultata  sulvants  : 

-  dans  la  figure  (9),  mesure  du  spectre  direct  par  la  perche  drolte  ;  IS  encore,  l'accord  avec  le  spectre 
de  Karman  est  excellent, 

-  dans  la  figure  (10),  mesure  de  la  coherence  entre  les  turbulences  mesurees  sur  les  alles  drolte  et 
gauche,  coherence  qul  d6cro!t  tres  nettement  en  fonction  de  la  frequence,  et  peraet  de  douter  qu'un 
controle  des  charges  en  boucle  ouverte  solt  efflcace  au-del4  de  0,6  Hz. 

Sur  la  me  me  figure,  on  a  aussl  trace,  un  polntllie,  la  coherence  pr£vue  th6orlquement  pour  une  tur¬ 
bulence  lsotrope  correspondant  au  spectre  de  Karman,  et  donn£e  par  la  fornule  : 


C(^)  = 


- ? -  (Md.\S/6  k  /  yd  \ 

r  (5/6)  V  2  V  '  a/6'  V  ‘ 


ou  r  est  la  fonction  eul£rienne,  d  la  distance  entre  les  mesures  de  turbulence,  V  la  vltesse  de  vol  et  K 
la  fonction  de  Bessel  modlfl£e  de  seconde  espSce. 

La  comparalson  "calcul-essal"  est  nettement  en  faveur  du  concept  d'lsotropie  de  la  turbulence. 

Ill  -  STATISTIQUES  DE  TURBULENCE 

III.l  -  Hlstorlque 

La  connalssce  statlstique  de  la  turbulence  atmosphSrique  a  6t6  obtenue  (et  est  tou jours  obtenue)  par 
le  blals  de  mesure  des  acc£l£ratlons  subles  par  les  avlons.  L' exploitation  de  ces  mesures  4  conduit  aux 
dispositions  rEglementaires  par  la  dSmarche  sulvante  : 

-  Etabllssement  d'une  mSthode  standard  de  calcul  de  la  rSponse  d'un  avion, 

-  Cholx  d'une  forme  de  rafale  standard, 

-  DSpoulllement,  gr4ce  4  ces  deux  modtles  des  accelerations  recuelllles,  c'est-4-dlre  obtentlon  des 
vlteases  de  rafales  en  valeur  et  en  nombre  par  distance  parcourue, 

-  Cholx  probablllste  de  la  vltesse  llmlte  de  rafale, 

-  Introduction  dans  la  reglementatlon  de  cette  valeur  llmite  et  du  mod41e  correspondant  de  calcul  de  la 
reponse  de  l’avion. 

Hlstorlquement,  cette  demarche  a  comportC  trols  etapes  : 

-  PCrlode  1931  4  1949  :  premieres  recherches  sur  une  rafale  lsoiee  tr4s  simple,  4  partlr  de  quelques 
r6sultats  de  vol  (travaux  de  Rhode  et  Donely), 

-  PCrlode  1950  4  1956  :  modlle  de  rafale  en  (1-coslnus)  et  travaux  de  Pratt,  Introduction  du  "Gust 
Alleviation  Factor"  ;  modile  d’ atmosphere  dCduit  de  vols  plus  nombreux  (55000  heures)  ;  etabllssement 
d'un  r4glement  qul,  pour  l'essentlel,  est  toujours  en  vigueur, 
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-  PErlode  1956  1  1970  :  travaux  de  Press  et  Houbolt  sur  la  turbulence  continue  reprEsentEe  par  un  spectre 
de  Kansan,  et  publication  par  la  FAA  en  1970  d'un  aaendeaent  dQ  2  la  synthiae  de  Hobllt. 

Cet  hlstorique  fait  clalreaent  apparattre  que  les  vltesses  de  rafales  qui  servent  de  base  2  la 
riglementation  ne  constituent  pas  une  description  de  la  rfialltE  de  l'atmosphEre,  nals  qu'elles  ne  sont  en 
fait  qu'un  moyen  de  reller  les  pica  d'accElEration  qu'on  peut  espErer  observer  sur  un  avion  donnE  aux 
pics  d'accfilEratlon  mesurEs  sur  un  autre  avion  pour  un  vol  dans  le  m£me  air  turbulent. 

Les  motivations  pour  une  reprise  d'activltE  dans  ce  domalne  ont  EtE  nombreuses  ;  nous  clterons  en 
partlculler  : 

-  la  nEcessltE  de  s'assurer  que  les  avlons  qui  volent  actuelleaent,  et  dont  on  cherche  2  augmenter  la 
durEe  de  vie,  sublssent  2  peu  pr2s  les  charges  pour  lesquellea  ils  ont  EtE  calculEs  ; 

-  le  besoln  de  valider  le  noddle  d' atmosphere  proposE  par  la  rEgleaentation  pour  un  beaucoup  plus  grand 
nombre  de  vols,  et  de  vErlfler  que  ce  noddle  n'est  pas  trop  sensible  au  type  d'avlon  sur  lequel  ont  EtE 
mesurEes  les  accElEratlons  ; 

-  1'lntErdc  d'une  extension  du  noddle  d'ataosphdre  aux  Evdnenents  trds  peu  probables  (trds  fortes 
rafales)  et  d'une  vErlflcatlon  que  ces  Evdnenents  aulvent  une  lol  de  Poisson  (probldme  de  la  rEslstance 
rfislduelle). 

Depuls  1980,  l'ONERA,  en  collaboration  avec  la  CAA,  a  abordd  ces  probldmes  par  la  crEatlon  d'une 
banque  de  donnEes  qui  porte  malntenant  sur  cinq  ans  d'activltE  de  la  flotte  de  British  Airways.  Air 
Prance  collaborera,  dds  1986,  2  l'Elarglssenent  des  informations  dlsponibles. 

111.2  -  La  banque  de  donnEes  de  British  Airways  : 

Depuls  quelques  annEes,  la  CAA,  par  le  progranae  CAADRAP ,  Etudlalt  systEmatlquenent  les  Evdnenents 
survenus  lors  de  vols  de  British  Airways,  pour  lesquels  un  facteur  de  charge  de  0,7  g  avalt  EtE  dEpassE. 
Ces  Evdnenents  Etant  trop  races  pour  l'Etabllssement  de  statlstlques  prEclses,  l'ONERA  a  demandE  que  le 
seull  solt  abalssE  2  0,5  g,  seull  pour  lequel  les  informations  seralent  en  nombre  ralsonnable,  nals  auf- 
flsant  pour  la  dEtermlnatlon  de  lols  statlstlques. 

L'ensemble  de  la  flotte  de  British  Airways  est  concern#  par  ce  collatlonneaent  de  donnEes, 
c'eat-2-dire  : 

9  Trident  2-E 
25  Trident  3-B 
21  BAC  1-1 l 
9  Trlstar  1011-1 
6  Trlstar  1011-200 
2  Tristar  1011-100 

16  Boeing  747-100 
11  Boeing  747-200 
6  Concorde 

28  Boeing  737 

17  Boeing  757. 

Au  cours  de  la  pErlode  de  5  ans  allant  du  ler  mal  1980  au  30  avril  1985,  environ  400000  vols  ont  EtE 
analysEs,  qui  correspondent  2  environ  800000  heures  de  vols  et  2  300000000  de  NM  parcourus  ; 
c'est-2-dlre  que,  compte  tenu  des  alEas,  60  X  des  vols  de  la  compagnle  ont  EtE  analysEs.  Au  cours  de 
cette  pErlode,  on  a  observE  environ  6300  Evdnenents,  dus  2  la  turbulence,  pour  lesquels  le  facteur  de 
charge  observE  dEpassalt  0,5  g. 

La  banque  de  donnEes  comprend  deux  parties  : 

a)  une  analyse  des  vols  par  avion  et  par  type  d'avlon  qui  donne  : 

-  le  nombre  d'heures  de  vol  effectuSes  chaque  mole  par  chaque  avion  d'un  type  donnE  ;  le  tableau  1 
prEsente  cette  information,  pour  les  50  premiers  mols,  pour  la  flotte  des  21  Bac  111, 

-  pour  l'ensemble  des  appareils  d'un  type  donnE,  et  par  mois,  le  nombre  d'heures  de  vol  et  le  nombre  de 
vols  analysEs,  le  nombre  de  vols  rEellement  effectuEs,  le  pourcentage  de  vols  analysEs  et  la  durEe 
moyenne  d'un  vol  (tableau  II), 

b)  une  llate  d' Informations  pour  chaque  Evdnement  observE,  qui  comprend,  par  colonnes  successlves 
(tableau  111)  : 

-  le  type  de  1 'avion, 

-  le  numEro  repEre  de  l'avlon  pour  le  type  correspondent , 

-  un  code  de  quatre  chlffres  donnant  les  conditions  de  vol, 

-  le  Jour,  le  mols,  1'annSe  de  I’EvEneaent, 

-  1'aEroport  de  dEpart, 

-  1'aEroport  d'arrivEe, 


-  la  aasse  de  l'avlon  (kg)  au  moment  de  1'lncldent, 

-  la  vltesse  lndlqu£e  (IAS), 

-  la  vltesse  vrale  (TAS), 

-  1 'altitude  (ft) 

-  la  position  des  volets, 

-  la  position  des  bees  de  bord  d'attaque, 

-  la  dur£e  de  l'dvdnement  (perlode  pendant  laquelle  des  accSlSratlons  super  leures  A  0,4  g  ont  ete 
observSes), 

-  le  nombre  d'heures  de  vol  de  la  cellule, 

-  le  nombre  de  Hach, 

-  le  C-  „  correspondent , 

-  le  facteur  de  charge  posltlf  (en  g), 

-  le  facteur  de  charge  ndgatif. 

III. 3  -  Restitution  de  la  turbulence  : 

La  banque  de  donnSes  creee  en  collabortlon  avec  la  CAA  a  ete  utilisSe,  de  dlffSrentes  maniSres,  en 
vue  d'obtenlr  des  statlstiques  de  turbulence  par  tranche  d'altltude.  On  s'est  alors  heurtS  au  probleme, 
Svoque  au  paragraphe  III. I,  de  trouver  un  aoyen  de  falre  correspondre  aux  pics  d’accSISratlon  des  pics  de 
turbulence.  L'approche  cholsle  a  6t6  celle  pr£conlsee  par  J.  Hall  (1962)  et  dScrite  par  J.  Taylor 
(pp.  208-210)  dans  son  "Manual  on  Aircraft  Loads"  ;  dans  cette  methode,  le  facteur  K  de  reponse  %  la 
turbulence  est  explicit®  en  nSgllgeant  le  mouvement  de  tangage  et  l'effet  de  1'envergure  flnle  sur  la 
portance  due  3  la  turbulence.  Un  programme  de  calcul  a  ete  realise  qul  permet  de  determiner  automatl- 
quement,  par  cette  methode  et  pour  chaque  type  d'avlon  la  turbulence  Squlvalente  et  la  turbulence  vrale 
(en  m/s). 


On  erfe  alors,  pour  chaque  type  d'avlon,  une  banque  de  donndes  partlculldre  dont  un  exemple  (relatlf 
au  Boeing  747)  est  donn®  dans  le  tableau  IV  ;  les  douze  premieres  colonnes  de  ce  tableau  sont  lea  semes 
que  celles  du  tableau  III  ;  on  a  ensulte  successlvement  : 

-  le  nombre  de  Mach, 

-  le  Cz  a  , 

-  le  facteur  de  charge  posltlf, 

-  le  facteur  de  charge  n®gatlf, 

-  la  turbulence  ®qulvalente  positive, 

-  la  turbulence  Squivalente  negative, 

-  la  turbulence  vrale  positive, 

-  la  turbulence  vrale  negative, 

-  la  dur®e  de  l'®v®nement, 

-  puls  les  num®ros  de  mols,  de  semaine  et  de  jour  depuls  le  d®but  du  collatlonnement  (ler  mal  1980). 

Partant  de  cette  banque  de  donn®e,  un  second  programme  calcule  alors  le  nombre  de  d®passements 
d'amplltudes  crolssantes  de  turbulence  pour  une  tranche  d'altltude  donn®e  ;  les  seuils  de  turbulence  3 
partlr  desquels  on  calcule  les  dSpassements  sont  de  6  m/s  3  16  m/s  pour  la  turbulence  Squivalente,  avec 
un  pas  de  0,3  m/s. 

Les  tranches  d'altldude  cholsles  pour  ces  comptages  sont  : 

o  <  z  <10000  ft 

10000  ft  <  z  <  20000  ft 
20000  n  <  z  <  30000  ft 
30000  ft  <  Z 

Si  on  connait  la  distance  parcourue  par  la  flotte  dans  la  tranche  d'altltude  consld®r®e,  on  en 
d®dult  ImmSdlatement  le  nombre  de  d®passements,  par  mile  nautlque  de  chaque  seuil  succesaif  de  comptage. 
C'est  cette  dernidre  information  qul  est  donnde  (sous  reserve)  dans  la  figure  11,  pour  l'ensemble  des 
Boeing  747  (types  100  et  200),  et  par  des  altitudes  aupdrieures  3  30000  ft  . 
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On  remarquera  que,  dans  la  representation  semi-log  choisle,  la  courbe  : 

E/NM  .  MW) 

est  une  droite  ;  c'est-3-dire  que  le  nombre  d'£v3nements  par  mile  nautlque  decrolt  exponentiellement  avec 
l'amplltude  de  turbulence,  r£sultat  coherent  avec  la  th€orle  des  valeurs  extremes. 


Ill .4  -  Analyse  critique  des  r6sultats  : 

Le  type  de  r6sultat  present!  dans  la  figure  11  a  £t!  calcul!  pour  chacune  des  families  d'avlons  utl- 
lls6es  par  British  Airways,  3  l'exception  des  Concorde  ;  dans  chaque  cas,  les  courbes  ont  senslblement  la 
mgme  allure  de  d€crolssance  exponentlelle  avec  l'amplltude. 

Par  contre,  si  on  reporte  sur  la  mgme  planche  le  nombre  de  d€passements  de  niveaux  de  turbulence  par 
mile  nautlque,  tel  qu'll  est  d£duit  des  accelerations  mesur£es  sur  dlff€rentes  families  d'avlons,  on  volt 
que  les  courbes  correspondantes  sont  loin  de  se  superposer  ;  c'est  ce  que  montre  la  figure  12  pour  la 
tranche  d* altitude  0-10000  ft. 

Deux  conclusions  peuvent  en  etre  tlr£es,  les  formules  simples,  telles  que  les  formules  de  Pratt  ou 
Hall,  sont  en  v£rit€  trop  simples,  et  1' information  qu'elles  fournlssent  sur  la  turbulence,  a  partlr  des 
accelerations  mesur£es  sur  une  famille  d'avlon  ne  permettent  pas,  3  l'lnverse,  de  prevoir  les 
accelerations  subies,  pour  la  meme  turbulence,  par  un  autre  type  d'apparells  ;  c'est  remettre  en  cause 
pratlquement  les  bases  de  la  certification  actuelle. 

Ces  formules  simples  sont  sufflsamment  correctes,  mals  certains  paramdtres  qu'elles  lmpliquent  ont 
£t£  mal  apprecies  (echelle  de  turbulence,  C^.a  Proflls  de  v°l)> 

PlutSt  que  de  remettre  en  cause  a  priori  l'utillsation  des  m£thodes  simples  qui  permettent  seules 
l'analyse  statlstlque  de  mllllers  d'6vdnements,  nous  nous  interesserons  lcl  aux  erreurs  syst£matlques 
susceptlbles  de  se  produlre. 

Nous  examinerons  successlvement  : 

a)  L' Influence  d'un  cholx  erron£  de  l'£chelle  de  turbulence  : 

Des  calculs  d'erreurs,  effectu£s  sur  la  formule  de  Hall,  ont  heureusement  montr!  une  faible  Influ¬ 
ence  de  ce  paramdtre. 

b)  L' Influence  d'une  mauvalse  appreciation  des  Cj,  a  : 


II  s'aglt  13  d'un  probldme  beaucoup  plus  grave,  car  toute  erreur  Bur  le  C^a  se  r£percute  entiS- 
rement  sur  le  calcul  de  la  turbulence,  et  les  Cz  ,  <X  sont  mal  connus  ;  on  dispose  en  effet  de  deux 
valeurs  du  Cz,a  fournies  par  les  constructeurs  :  l'une,  le  Cz.GL  dit  "rlglde"  provlent  directement  des 
mesures  en  soufflerle,  tandis  que  l'autre,  le  Cz,a  dlt  “flexible"  en  est  d£dulte  par  un  calcul  qul  prend 
en  compte  les  deformations  aeroeiastlques  quasl-statlques  de  la  structure.  Le  cholx  entre  les  deux  depend 
des  hypothdses  que  l'on  est  amen!  3  faire  sur  la  bri3vet£  de  la  rafale. 

La  figure  13  pr£sente  les  nombres  de  dCpassements  de  la  turbulence  calcul€s  3  partlr  des  accele¬ 
rations  mesur£es  sur  les  Boeing  747  dans  la  tranche  d'altltude  0-10000  Ft  ,  sulvant  les  deux  hypotheses 
sur  les  Cz  a*  Lm  differences  d'allure  sont  manlfestes. 

c)  Les  proflls  de  vol  : 

Les  seules  Informations  fournles  par  la  CAA  sur  l'utillsation  des  avlons  sont  les  nombres  d'heures 
de  vol  et  le  nombre  de  vols  analyses  chaque  mols.  On  Ignore  done  a  priori  le  nombre  d'heures  analys£es 

(oO  le  nombre  de  miles  nautlques  parcourus)  pour  chaque  tranche  d'altltude.  Une  analyse  statlstlque  fine 
des  missions  de  la  flotte  et  l'examen  detain!  d'un  certain  nombre  de  vols  repr£sentatlfs  ont  permls  de 
dCflnlr,  pour  les  Boeing  747,  un  profll  de  vol  moyen  tr3s  probablement  proche  de  la  r£allt£.  Pour  les 
autres  families  d'avlons,  on  dispose  seulement  d' informations  f ragmentalres,  et  les  proflls  de  vol 
cholsls  pour  l'analyse  peuvent  tr3s  blen  s'£v£rer  faux  dans  un  rapport  de  1  3  2  pour  certalnes  tranches 
d'altltude. 

Ill .5  -  Distribution  statlstlque  des  paquets  de  turbulence  : 

Alnsl  qu'll  a  £t£  dlt  dans  1' Introduction,  on  s'est  lnt£ress£  3  la  distribution  statlstlque  des 
paquets  de  turbulence,  distribution  lnt£ressante,  tant  du  point  de  vue  d'une  mod£llsatlon  de  l'atmosphdrc 
que  de  celul  de  probldmes  lies  3  la  residence  rCslduelle.  On  est  parti  de  la  remarque  que,  si  les  rafales 
Ctalent  surement  corr£l£es  entre  elles  dans  un  paquet  de  turbulence  (existence  de  noddies  spectraux),  11 
n'y  avalt  aucune  raison  que  les  paquets  de  turbulence  eux-mgmes  ne  solent  pas  en  g£n£ral  Independants. 
Dans  ces  conditions,  la  distribution  de  paquets  de  turbulence  o3  une  rafale  svalt  d£pass£  un  seull 
devalt  obClr  3  une  lol  de  Poisson. 

Pour  vCrlfler  cette  hypothise,  et  prenant  pour  unit!  de  temps  la  semalne,  on  a  proc£d£  de  la  manldre 
sulvante  (l'exemple  du  tableau  V  tralte  de  turbulences,  entre  0  et  10000  ft  d£dultes  de  vols  de  la 
flotte  des  Boeing  747). 

On  a  compt£  le  nombre  de  paquets  de  turbulence  pour  lesquels  un  niveau  de  11  m/s  avalt  £t£  dispose 
(137  pour  cet  example)  ;  on  connalt  par  allleurs  le  nombre  N  de  semalnes  soumlsea  3  analyse  (lcl 
tt  "  241*  ;  on  en  dCdult  le  nombre  moyen  |L  d'Cvdnements  par  semalne.  Appllquant  la  formule  assoclCe  3  la 
distribution  de  Poisson  : 
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Prob  (K  SvSneoents) 


on  connait  alors  la  probability  pour  qu'une  semaine  appartlenne  3  la  classe  0  (0  yvdnement),  3  la 
clasae  1  (1  yvSnement),  3  la  classe  K  (K  €v3nements) . . . .  Le  programme  Sexploitation  compte  par 

allleurs  le  nombre  N(K)  de  semalnes  appartenant  3  la  classe  K.  On  en  dSdult  la  probability 
"expyrlmentale"  : 


Prob  (K  SvSnements) 


N(l<) 

N 


Le  tableau  V,  qul  prysente  la  comparaison  entre  la  probability  thyorique  et  la  probability  obaervSe 
pcur  chaque  classe  montre  3  l’evldence  que  la  distribution  de  Poisson  est  reprysentatlve  de  la  ryality. 

IV  -  CONCLUSIONS 

Un  tr3s  gros  effort  a  ?te  fait  au  cours  de  ces  dernieres  annees ,  tant  dans  le  domalne  des  mesures 
fines  de  la  turbulence  que  dans  celul  de  l'analyse  statlBtique  des  masses  d' Informations  obtenues  sur 
avion*  commerciaux  au  cours  de  vols  ryguliers.  Les  conclusions  de  ce  rapport  sont  nettes  sur  le  premier 
point,  mals  conslderablement  plus  nuancyes  sur  le  second. 

a)  Les  mesures  fines  de  turbulence  sur  Mirage  3  et  Nord  260  spyclalement  Iqulpes  ont  montre  que  l'on 
avalt  acquit  la  maltrlse  de  la  mesure  en  temps  ryel  de  rafales  vertlcales  dans  une  grande  plage  de 
fryquences  (plus  de  12  Hz).  Les  turbulences  mesurles  ont  des  densltys  spectrales  qul  correspondent 
parfaitement  au  module  de  Karman  ;  les  cohyrences  entre  turbulences  mesuryes  en  dlffirents  points  confor- 
tent  par  allleurs  le  concept  d'lsotrople,  ou  tout  au  molns  d'orthotrople  dans  le  plan  horizontal.  On 
remarque  enfin  que  les  methodes  de  calcul  des  ryponses  de  l'avlon  aux  rafales  ont  permis  une  bonne  pryvi- 
sion  des  rysultats  d'essals  en  vol,  meme  quand  les  modes  souples  etalent  concernys. 

L' effort  dolt  porter  malntenant  sur  la  mesure  en  vol  de  la  turbulence  transverse,  opyration  beaucoup 
plus  dSllcate. 

b)  Grace  3  la  CAA  et  3  British  Airways,  l'ONERA  a  pu  cier  la  plus  grande  banque  exlstante  (3  notre 
connalssance)  sur  les  facteurs  de  charges  dus  3  la  turbulence  ;  le  choix  d'un  seull  de  0,5  g  rysulte 
d'un  compromls  entre  la  finesse  d'informatlon  souhalt?le  et  les  possibilites  d'exploitation.  Le  collation- 
nement  des  donnyes  sera  poursulvl,  et  Air  France  ae  joindra  au  club  en  Janvier  1986. 

L'exploitation  des  donnyes  prysente  de  syrleuses  difficulty*  si  les  depassements  par  tranche  d'altl- 
tude  d'un  seull  donny  de  turbulence  sont  cohyrentes  3  l'intyrieur  d'une  famllle  d'avlons,  il  n'en  est  pas 
de  mSme  quand  on  r6unlt  les  informations  dydultes  des  diffyrentes  families.  Cette  difficulty  peut  Stre 
liye,  solt  au  fait  que  l'utlllsatlon  de  formules  simplifies  de  calcul  de  la  turbulence  n'a  pas  de  sens 
(la  philosophle  de  la  certification  seralt  alors  ybranie),  soit  au  fait  que  les  paramdtres  lntroduits 
dans  ces  formules  sont  mal  appryciys,  II  est  pour  l'instant  impossible  de  trancher  entre  ces  deux  hypo¬ 
theses. 

On  remarquera  enfin  que  les  paquets  de  turbulence  obyissent,  co me  pryvu,  a  une  distribution  de 
Poisson. 

Nous  pensons  que  1'effort  dolt  porter  partlcull3rement  sur  les  points  sulvants  : 

-  Evaluation  de  nouvelles  formules  slmpllfiyes  de  restitution  de  la  turbulence  rycemment  proposyes 
(la  figure  16  donne  une  comparaison  entre  les  dypassements  de  turbulence  dydults  des  mymes  acc6 irations 
par  la  formule  de  Pratt  et  la  formule  de  Turner). 

-  Amelioration  de  la  connalssance  des 

-  A  cours  d'essals  en  vol,  mesures  slmultanyes  de  l'accyiyration  et  de  la  turbulence  dont  on  dydulra 
le  rapport  entre  pics  d'accyiyratlon  et  pics  de  turbulence  ;  comparaison  de  ce  rysultat  experimental  avec 
celul  pryvu  par  diffyrentes  formules  slmpllfiyes. 

Les  probldmes  abordys  lcl  sont  lmportants  sur  plus  d'un  aspect,  et  concernent  l'ensemble  de  la 
communauty  atlantlque.  Une  mlse  en  commun  des  donnyes,  une  comparaison  des  rysultats  et  une  critique 
ryclproque  semblent  indlspensables. 
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ABSTRACT 

There  has  been  considerable  interest  in  the  U.K.  over  a  long  period  in  the  response  of  an  aircraft  in  turbulence,  and  in 
the  analysis  of  turbulence  itself.  Work  is  discussed  which  includes  the  collection  and  analysis  of  operational  flight  data, 
methods  used  to  calculate  aircraft  response  in  turbulence,  including  non-linear  response,  and  low  level  operation.  This  paper 
provides  an  overview  of  the  U.K.  work  in  this  field. 


1.0  INTRODUCTION 

When  it  was  discussed  in  San  Antonio  in  April,  1 985,  at  an  ad-hoc  meeting  of  those  present  at  the  Structures  and 
Material  Panel  meeting  who  were  interested  in  turbulence,  it  was  suggested  that  there  were  three  topics  of  particular 
interest:— 

(i)  Collection  and  analysis  of  flight  data 

(ii)  Methods  of  calculating  aircraft  response  in  turbulence  (including  non-linear  response) 

(iii)  Low  level  operation  (including  correlation  of  ground  contours  with  turbulence). 

Following  the  San  Antonio  meeting,  discussions  in  the  U.K.  led  to  an  informal  meeting  where  the  work  in  the  U.K.,  both 
in  government  establishments  and  industry,  was  reviewed.  This  work  is  presented  briefly  in  this  paper  which  thus  provides  an 
overview  of  the  U.K.  work  and  serves  to  present  the  various  topics  on  which  work  has  been  done  and  is  currently  underway, 
and  the  topics  on  which  work  is  planned  in  the  future.  Material  is  available  on  all  three  subjects  noted  above  and  it  is  possible 
to  identify  a  number  of  papers  which  could  be  presented  at  a  future  Specialist  Meeting  on  Turbulence. 

The  following  paragraphs  have  in  general  been  penned  by  those  active  in  the  various  topics,  and  acknowledgement  of 
those  responsible  is  made  at  the  end  of  the  paper. 


2.0  COLLECTION  AND  ANALYSIS  OF  OPERATIONAL  FLIGHT  DATA 

2. 1  Reassessment  of  Gust  Statistics  in  the  Light  of  Current  Analysis  Methods  (Propellor  Driven  Aircraft) 

Current  civil  gust  load  requirements  are  largely  based  upon  statistical  data  collected  from  U.S.  domestic  operations 
prior  to  1950  and  supported  by  data  measured  during  U.K.  and  European  operations  prior  to  1 960.  Gust  velocity  data  was 
derived  from  exceedance  counts  of  c.g.  acceleration  using  a  simple  gust  alleviation  factor  approach  based  upon  average  wing 
loading  data  and  assuming  a  rigid  aircraft.  This  approach  necessarily  employed  many  simplifying  assumptions  which  have 
been  eliminated  from  modem  gust  response  calculations.  The  gust  velocities  derived  from  this  accumulated  acceleration  data 
have  therefore  been  reviewed  in  the  light  of  current  analytical  techniques.  The  following  factors  have  been  investigated  :- 

(i)  The  effects  of  aircraft  flexibility  on  the  dynamic  response  characteristics  of  the  data  collecting  aircraft. 

(ii)  The  effects  of  contributions  from  fuselage  and  tailplane  on  the  aircraft  normal  force  coefficient. 

(iii)  The  effects  of  the  pitch  degree  of  freedom  on  the  dynamic  response  of  the  data  collecting  aircraft. 

(iv)  The  effects  of  allowing  for  time  delay  between  the  various  lifting  surfaces  penetrating  the  gust  front. 

(v)  The  sensitivity  of  the  response  of  the  data  collecting  aircraft  to  unsteady  lift  assumptions. 

The  investigations  have  re-evaluated  101,789  hours  of  acceleration  data  collected  on  three  different  types  of  transport 
aircraft  between  1 933- 1 958.  (References  2.1.1  and  2. 1 .2). 
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Results  show  that  many  of  the  refinements  which  have  been  introduced  into  the  mathematical  modelling  methods  have 
built  additional  conservatism  into  the  derived  gust  loads.  Thus  the  simplified  models  used  previously  to  assess  operational 
gust  statistics  have  led  to  conservative  estimates  of  derived  gust  exceedances.  It  is  also  shown  that  current  aircraft  modelling 
techniques  have  introduced  additional  increases  in  predicted  aircraft  response  loads  without  allowances  for  compensating 
reductions  in  design  gust  velocities. 


2.2  Further  Reassessment  of  Gust  Statistics  using  CAADRP  Data 

A  further  reassessment  of  Gust  Statistics  using  CAADRP  (Civil  Aviation  Authority  Data  Reduction  Programme)  Data 
has  been  made  by  analysing  the  response  of  a  BAC  1  - 1 1  500  'Series  aircraft'  on  standard  airline  routes. 

The  assumptions  and  procedures  used  in  producing  a  mathematical  model  for  the  1-11  500  series  aircraft  are  described 
and  an  assessment  of  the  measured  incremental  centre  of  gravity  (C.G.)  accelerations  from  CAADRP  and  ONERA  data  is 
made.  The  CAADRP  results  refer  to  a  set  of  routine  airline  route  results  using  all  flight  data  whereas  the  ONERA  data 
consisted  only  of  a  set  of  high  ‘g’  CAADRP  results  produced  for  incremental  accelerations  in  excess  of  0.5g. 

An  assessment  has  been  made  of  the  ‘average  mission"  from  each  of  the  1 24  routine  C  AADRP  flights  and  this  has  been 
broken  down  for  all  the  flights,  for  flights  in  the  Berlin  Corridor  and  for  flights  not  in  the  Berlin  Corridor.  Mission  Analysis 
results,  using  the  characteristic  values  of  frequency  (N„)  and  Root  Mean  Square  (A),  are  determined  from  the  mathematical 
model  of  the  1-11  500  series  aircraft  calculated  in  accordance  with  current  airworthiness  regulations,  with  the  measured 
CAADRP  and  ONERA  data  and  these  have  been  presented  in  reference  2.2.1.  Also  included  in  the  Appendix  is  a 
comparison  of  the  Mission  Analysis  results  with  those  using  the  ESDU  "Average  Gust  Frequencies'  method  presented  in 
reference  2.2.2.  This  comparison  has  been  made  for  incremental  C.G.  accelerations.  Comparisons  have  also  been  made  for 
all  mission  segments  and  the  average  design  cruise  speed  comparing  the  incremental  C.G.  accelerations  obtained  from  the 
Mission  Analysis  method  with  those  of  Pratt,  discrete  1  -cosine  gust  and  design  envelope  methods. 

From  the  results  presented  for  the  above  methods  conclusions  have  been  drawn  comparing  predicted  gust  loads  for  a 
probability  of  2x  1 0“'  exceedances  per  hour  (as  defined  in  JAR  25.305(d))  with  those  found  from  the  measurements 
contained  in  the  CAADRP  and  ONERA  data  for  a  1-11  500  series  aircraft. 

The  gust  velocities  obtained  from  the  1-11  ONERA  data  have  also  been  compared  with  the  values  obtained  from  the 
original  assessment  of  the  American  data  based  on  measurements  prior  to  1 950  and  on  which  certification  requirements  for 
gust  loads  are  based. 


2.3  Gust  Description  for  Active  Control  Technology 

The  design  of  aircraft  using  active  control  technology  requires  a  knowledge  of  the  individual  gusts  that  make  up 
continuous  turbulence.  A  method  of  describing  the  shapes  of  the  individual  gusts  has  been  developed  from  data  produced  on 
research  flights  on  a  T33  at  1000ft  and  on  a  Canberra  near  the  top  of  a  thunderstorm  at  46000ft.  It  is  shown  that  the  gusts 
that  are  important  for  a  particular  response  are  strongly  dependent  on  the  rate  at  which  Zero  Crossings  occur  in  that 
response.  The  method  has  been  used  on  CAADRP  special  turbulence  events  measured  on  the  BAC  1-11  and  the  Trident. 
(Reference  2.3.2). 

A  plot  has  been  made  of  peak  magnitude  versus  wave  length.  This  is  done  by  plotting  in  a  non-dimensional  form  peak 
magnitude  divided  by  standard  deviation  and  wave  length  divided  by  average  wave  length  1  /N,„  where  N„  is  the  count  of 
zero  crossings. 

The  shape  of  this  curve  is  shown  to  be  very  similar  whether  it  be  from  gust  velocity  data  or  C.G.  data  from  T33, 
Canberra  or  CAADRP  record,  and  whether  from  full  or  reduced  sampling  rates. 

Finally  a  curve  can  be  drawn  for  all  the  CA  ADRAP  data  on  which  is  drawn  a  curve  deduced  from  the  points  at  a 
probability  of  1  in  1 00  augmented  by  an  allowance  for  manoeuvring  equivalent  to  1  m/ sec. 

This  curve  may  be  said  to  be  a  result  that  the  work  was  seeking,  representing  the  gust  magnitude  as  a  function  of  gust 
wave  length  for  a  constant  global  probability.  The  global  probability  is  deducible  from  the  1  in  1 00  of  the  figure  itself 
multiplied  by  the  probability  of  a  CAADRAP  Special  Event  of  this  severity  assessed  at  about  once  every  1 1 ,000  to  1 8.000 
hours  which  implies  limit  load  levels. 

It  is  concluded  (in  Reference  2.3.3)  that  the  results  presented  should  suffice  to  enable  a  relaxation  of  the  isolated  gust 
regulation  at  the  lower  wave  lengths. 

Further  work  is  planned  looking  into  the  significance  of  combined  loads  from  gusts  and  manoeuvres,  and  the  control 
movements  which  should  be  allowed. 


2.4  The  Operating  Environment  for  Modern  Transport  Helicopters 

The  CAA.  in  concert  with  major  rotorcraft  operators  and  ha»  'w"‘"nllv  laiiflfitoi  *  dill  ttBaafaft 


programme  on  two  types  of  transport  helicopter.  It  is  intended  that  strain  gauge  and  accelerometer  data  measured  during 
routine  civil  operations  will  provide  a  more  detailed  knowledge  of  the  service  operating  environment  for  these  unique  types 
of  aircraft.  Of  particular  interest  will  be  the  collection  of  additional,  low  level  gust  data  for  comparison  with  existing  statistics 
derived  entirely  from  fixed  wing  operations.  It  is  hoped  that  sufficient  data  will  be  available  by  the  end  of  1 986  to  allow 
examination  of  the  results  and  an  assessment  of  the  implications  on  the  response  of  fixed  wing  aircraft. 

2.5  Measurement  of  the  Fore  and  Aft  Component  of  Turbulence  on  Operational  Combat  Aircraft 

An  instrument  to  record  the  high  frequency  fluctuations  in  pitot  pressure  caused  by  fore  and  aft  gusts  experienced  by 
aircraft  has  been  designed.  It  is  being  fitted  to  a  few  combat  aircraft  which  will  begin  collecting  data  in  RAF  operational 
service  in  the  Autumn  of  1985.  The  aircraft  are  also  instrumented  to  record  loads  experienced.  Techniques  for  handling  and 
reducing  the  turbulence  data  are  at  present  being  evolved  and  it  is  hoped  that  a  statement  on  progress  can  be  made  by  the 
end  of  1986. 


3.0  METHODS  OF  CALCULATING  AIRCRAFT  RESPONSE  IN  TURBULENCE  (INCLUDING  NON-LINEAR 
RESPONSE) 

3. 1  Modelling  the  Flexible  Aircraft  and  the  Turbulence 

A  brief  description  is  given  of  the  way  in  which  the  flexible  aircraft  and  the  turbulence  are  modelled.  The  linear 
representation  which  follows  has  been  in  existence  for  many  years.  The  non-linear  work  is  of  recent  origin. 

The  basic  aircraft  equations  of  motion  are  assembled  as  a  set  of  simultaneous  second  order  differential  equations.  The 
generalised  co-ordinates  are  chosen  as  normal  modes  or  branch  modes  of  the  aircraft,  and  the  response  forces  are  equated  to 
the  gust  forces.  The  mode  shapes  are  chosen  to  cover  the  frequency  range  of  interest  and  the  response  quantities  required. 

The  equation  of  motion  may  be  solved  for  a  specific  gust  input  by  either  of  two  principal  methods.  The  time  plane 
solution  which  involves  a  step-by-step  integration  of  the  equation  timewise  for  the  required  input,  and  the  frequency  plane 
solution  which  involves  derivation  of  a  transfer  function  for  the  response  quantity  of  interest  to  a  unit  gust  input.  This  transfer 
function  is  then  multiplied  by  the  input  gust  spectrum,  and  an  output  response  spectrum  obtained,  from  which  may  be 
calculated  the  continuous  turbulence  co-efficients  familiarly  A  and  N„. 

Both  forms  of  solution  are  appropriate  where  the  representations  of  the  aircraft  and  aerodynamics  are  linear,  and 
Fourier  transform  methods  may  be  used  to  show  the  exact  equivalence  of  the  two  approaches.  The  standard  calculations 
carried  out  for  Airworthiness  and  Procurement  Authorities  assume  a  standard  linear  representation,  and  Continuous 
Turbulence  (C.T.)  requirements  have  been  formulated  on  a  statistical  basis  using  linear  equations. 

3.2  Response  Solutions  for  Non-Linear  Representations 

Non-linearities  can  be  introduced  in  either  the  aircraft  modelling  or  the  aerodynamics,  or  both.  Non-linearities 
normally  encountered  include  non-linear  lift  curve  slope,  non-linear  power  control  units  or  control  surfaces,  digital  control 
and  non-linear  control  laws. 

Whereas  the  time  plane  solution  may  be  used  with  no  inherent  difficulty  in  the  non-linear  case,  the  assumptions 
required  for  the  frequency  plane  method  are  immediately  violated. 

Solutions  for  discrete  gust  inputs  present  no  problems  using  the  time  plane,  but  solutions  for  C.T.  have  to  cover  long 
periods  of  time  to  give  statistically  reliable  results.  Work  is  in  hand  to  develop  acceptable  approximate  methods  of  producing 
C.T.  results. 

The  initial  approach  is  to  produce  pseudo-random  gust  histories  against  time  which  when  applied  to  a  linear  time  plane 
aircraft  model  will  give  results  in  close  agreement  with  results  obtained  for  a  conventional  frequency  plane  C.T.  analysis. 
These  gust  histories  may  then  be  used  as  the  input  to  the  non-linear  model,  and  response  quantities  produced.  In  practice, 
care  is  taken  to  minimise  the  effects  of  transients  due  to  initial  conditions,  and  to  obtain  a  satisfactory  resolution  at  low 
frequency  without  computing  for  an  undue  length  of  time. 

This  approach  can  be  made  to  work  but,  even  with  approximations,  is  extremely  cumbersome  to  use  and  excessively 
expensive  in  computer  time.  Even  using  a  reduced  model  only  a  limited  range  of  response  quantities  and  flight/aircraft 
conditions  can  be  covered  in  a  reasonable  time. 

To  cover  the  complete  range  of  cases  necessary  for  an  aircraft  design,  alternative  methods  are  being  tried  whereby 
discrete  gust  results  are  used  to  estimate  the  C.T.  results.  Early  results  look  promising  although  in  the  long  run.  these  results 
will  have  to  be  examined  together  with  the  results  obtained  using  the  statistical  discrete  gust  method,  and  their  adequacy 
judged. 
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3.3  Test  Techniques  for  Validation  of  Aircraft  Response 

There  is  a  need  to  validate  the  dynamic  and  static  mathematical  models  of  the  structure,  the  aerodynamics,  the  systems 
and  the  input  gust  environment  and  as  a  result  test  techniques  have  to  be  derived  that  cover  each  model  in  turn  and  then  the 
total  aircraft  and  system  in  a  flight  test. 

Testing  includes  ground  resonance  testing,  ground  rig  tests,  simulators,  wind  tunnel  tests  and  flight  tests.  Flight  testing 
includes  both  defined  control  tests  to  excite  the  aircraft  in  steady  conditions  and  measurements  of  the  response  in  turbulence. 

Recently  the  A3 10  and  the  BAC  1-11  have  been  used  for  flight  testing  load  alleviation  systems  and  a  short  summary  of 
the  test  techniques  are  given  below  :- 

•  A3 1 0  Load  Alleviation  Flight  Testing 

In  this  test  spoilers  and  outboard  ailerons  were  used  to  provide  the  alleviation,  and  the  testing  consisted  of  measuring 
the  transfer  functions  between  the  controls  and  a  number  of  interesting  quantities  for  the  purpose  of : 

—  checking  the  prediction  processes  of  the  transfer  function 

—  checking  the  stability  of  the  system  and  comparison  with  predictions 

—  checking  the  control  law  tolerances  in  relation  to  the  uncertainities  in  the  prediction  process 

—  checking  the  open  loop  time  histories  response  to  direct  ramp,  triangular  and  pseudo  continuous  turbulence  in  the  open 
and  closed  loop  of  the  load  alleviation  law  against  predictions 

These  tests  were  carried  out  at  two  speeds  and  the  tests  showed  that  the  functioning  of  the  load  alleviation  law  was 
successfully  demonstrated  in  flight  and  that  the  response  of  the  major  modes  of  the  aircraft  to  sinusoidal  and  discrete  inputs 
compared  reasonably  well  with  the  predictions. 

•  BAC  1-1 1  —  Load  Alleviation  Flight  Testing 

In  this  test  only  the  spoiler  controls,  using  fast  rate  jacks  (300*/sec),  were  used  to  provide  the  alleviation  for  the  wing. 

In  the  open  loop  a  large  number  of  flight  tests  were  carried  out  at  a  number  of  speeds  in  which  discrete  and  continuous  inputs 
were  put  into  the  control  system  over  a  wide  range  of  wavelengths,  amplitudes  and  frequency.  The  response  of  the  aircraft 
was  then  analysed  in  order  to  check  the  spoiler  aerodynamic  forces,  lag,  and  responses  of  the  aircraft  by  way  of  accelerations 
and  forces. 


3.4  Response  of  Flexible  Wind  Tunnel  Models  to  Turbulence 

The  U.K.  is  participating  in  the  design  and  testing  of  active  control  systems  for  flutter  suppression  and  gust  load 
alleviation  on  a  wind  tunnel  model.  The  model  is  basically  a  Tornado  flutter  model  which  has  been  modified  for  these  tests 
by  MBB  to  include  inboard  and  outboard  control  surfaces  on  the  trailing  edge  of  the  wings,  and  now  spoilers  on  the  upper 
surface  of  the  wing. 

Specific  turbulence  forms  have  been  generated  in  the  wind  tunnels  by  the  use  of  moving  vanes  and  the  responses 
measured.  In  comparison  with  flight  measurements,  the  lack  of  realism  in  some  aspects  of  the  experiment  is  offset  by  the 
unique  ability  to  generate  repeatable  discrete  gusts  and  continuous  turbulence  with  known  characteristics.  Aerodynamic 
representation  of  the  aircraft  and  the  effects  of  turbulence  penetration  were  investigated  and  the  mathematical  model  of  the 
aircraft  was  refined.  On  this  basis  an  active  control  system  was  designed  which  successfully  reduced  the  wind  tunnel  gust 
loads,  see  reference  3.4.1. 

In  the  U.K.  preliminary  tests  have  been  carried  out  in  the  British  Aerospace,  Brough  high  speed  wind  tunnel  where  a 
paddle  wheel  has  been  installed  in  the  tunnel  floor  to  act  as  a  gust  generator.  To  date  sinusoidal  testing  has  been  carried  out. 
and  work  is  planned  on  step  and  ramp  inputs,  see  reference  3.4.2. 


3.5  Spanwise  Variation  of  Turbulence 

Wing  loads  predicted  from  continuous  turbulence  models  are  compared  for  power  spectra  assuming  the  same  gust 
velocity  across  the  span  of  the  aircraft  and  spectra  including  spanwise  variations  of  gust  velocity.  Investigations  are  extended 
to  aircraft  with  active  controls  for  gust  load  alleviation.  It  is  shown  that  there  are  significant  differences  in  the  predicted 
alleviation  in  some  circumstances.  This  work  is  still  underway  at  RAE  Famborough  and  will  be  the  subject  of  a  report  still  to 
be  issued. 

Measurements  of  the  spanwise  variation  of  the  normal  component  of  turbulence  are  available  from  GNAT  XP505. 
Some  records  have  been  analysed  using  both  Power  Spectral  Density  and  Statistical  Discrete  Gust  techniques,  and  further 
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investigations  of  the  data  are  planned.  (Reference  3.5) 


3.6  SDG/PSD  Relationships  (Linear  Response) 

The  power-spectral-density  (PSD)  and  statistical-discrete-gust  (SDG)  methods  have  in  the  past  been  contrasted  as 
alternative  approaches  to  the  analysis  and  prediction  of  aircraft  response  to  atmospheric  gusts  and  turbulence.  Recent  work 
has  clarified  the  relationship  between  the  methods  and  provided  a  common  framework,  in  the  form  of  a  statistical  pattern 
theory,  whereby  a  unified  approach  (3.6.1)  may  be  adopted. 

In  a  recent  series  of  memoranda  (3.6.2-4),  it  has  been  shown  how  standard  power-spectral  procedures  to  quantify  the 
stochastic  response  of  linear  systems  may  be  implemented  by  a  variational  technique :  the  ‘method  of  equivalent 
deterministic  variables'.  In  Reference  3.6.3  a  discrete-element  technique  is  described  for  implementing  numerically  the 
above  variational  formulation.  The  input  profiles  used  in  this  analysis  are  already  familiar  from  existing  discrete-gust 
methods.  On  this  basis  it  has  been  shown  (3.6.3)  that  the  SDG  and  PSD  methods  may  be  matched  such  that  the  former  is 
essentially  an  approximate  numerical  implementation  of  the  latter.  On  the  other  hand,  reference  3.6. 1  also  discusses 
important  aspects  in  which  the  SDG  and  PSD  models  differ.  This  divergence  applies  particularly  to  the  modelling  of  gusts  of 
extreme  intensity. 


3.7  SDG  and  Non-Linear  Response 

A  feature  of  the  SDG  method  is  that  non-linear  and  linear  systems  may  be  assessed  in  a  unified  manner  3.7.1, 3.7.2.  In 
each  case,  response  evaluation  is  based  upon  a  worst-case  analysis  in  which  the  maximum  response  to  families  of  equi- 
probable  gust  patterns  is  found.  For  non-linear  systems,  however,  the  search  for  a  worst-case  must  be  performed  for  several 
families  of  inputs,  corresponding  to  different  levels  of  intensity,  whereas  for  linear  systems  only  a  single  (arbitrary)  intensity 
need  be  considered. 

The  validity  of  the  SDG  method  for  predicting  the  stochastic  response  of  non-linear  systems  has  been  supported  by 
numerical  simulation  studies  (3.7.3)  in  which  the  results  of  the  prediction  method  have  been  compared  with  measurements  of 
the  simulated  stochastic  response  of  various  modelled  systems.  In  this  digital-simulation  study,  inputs  were  random-process 
samples  and  system  outputs  were  calculated  by  numerical  integration  of  the  associated  non-linear  differential  equations. 

Applications  of  the  SDG  method  to  investigate  effects  of  non-linear  control-system  behaviour  upon  aircraft  response 
are  reviewed  in  Reference  3.7.2. 


3.8  Correlation  of  Stochastic  Loads 

A  particular  problem  arising  in  the  specification  of  structural  design-load  conditions  concerns  the  relationship  between 
the  results  of  a  PSD  statistical  analysis  of  the  aircraft  dynamic  response  and  the  requirements  dictated  by  existing  Stress 
Office  practice.  Stress  analysis  utilises  design  conditions  incorporating  a  set  of  forces  in  equilibrium.  PSD  methods  however, 
lead  instead  to  individual  design-level  values  of  load  such  as  shear,  bending  moment  and  torsion,  at  various  points  in  the 
structure.  They  are  defined  statistically  and  generally  do  not  occur  simultaneously. 

Various  methods  (e.g.  Reference  3.8.1)  exist  for  generating  design  conditions  to  match  the  statistically  defined  loads 
resulting  from  the  PSD  analysis.  More  recently,  a  rational  method  for  generating  design  conditions  which  incorporate  forces 
in  equilibrium  has  been  presented  by  Noback  (3.8.2).  This  method,  the  equal  probability  technique,  is  based  on  the  idea  that 
the  correct  design  load  condition  for  any  load  or  stress  is  that  combination  of  loads  that  is  most  probable  under  the  condition 
that  the  specified  load  or  stress  is  equal  to  its  design  value. 

Whilst  Noback’s  procedure  is  theoretically  ideal,  evaluation  of  the  large  number  of  correlations  involved  could  be 
prohibitive.  In  reference  3.8.3  an  approximate  method  was  therefore  proposed  which  combines  the  objectives  of  the  equal 
probability  technique  with  a  representation  of  design  loads  as  a  sum  of  elementary  load  distributions.  An  associated 
matching  process  was  also  proposed.  The  method  involves  the  use  of  deterministic  input  patterns  comprising  'ramp-gust' 
components  and  is  based  on  relationships,  recently  established,  (reference  3.6.1),  between  the  SDG  and  PSD  models  of 
turbulence  (see  para.  3.6). 


4.0  LOW  LEVEL  OPERATION 

4.1  Measurement  and  Analysis  of  Turbulence  using  an  MoD  GNAT  at  Bedford 

Low  altitude  turbulence  measurements  made  by  a  Gnat  aircraft  have  been  described  in  reference  4.1.1.  Two  types  of 
statistical  analysis  of  the  turbulence  records  have  been  employed :  power  spectral  density  and  statistical  discrete  gust.  In 
illustration  of  these  methods,  six  samples  of  turbulence,  all  at  mean  radio  altitudes  at  or  below  200  metres,  have  been 
discussed  in  detail  (4.1.1).  Analysis  of  a  much  larger  sample  of  the  turbulence  data  gathered  is  continuing  and  results  are  to 
be  published  in  a  later  report. 


A  particular  outcome  of  the  analysis  is  a  two-parameter  representation  (4.1.2)  of  turbulence  patches  which  combines 
measures  of  power-spectral-density  and  intermittency  with  an  exponential  (non-Gaussian)  amplitude  distribution.  This 
representation  is  applicable,  in  particular,  to  the  prediction  of  loads,  ride-bumpiness  and  control-system  activity. 

Another  particular  area  of  analysis  has  been  the  comparison  of  statistics  at  the  most  extreme  intensities  with  those  at 
lower  levels  of  intensity  (4. 1.3).  The  results  tend  to  support  a  recent  proposal  (4. 1 .4)  for  a  turbulence  model  in  which  the 
fractal  dimension  (4.1.S)  of  the  turbulence  takes  different  values  under  conditions  of  extreme  and  lesser  intensity. 


4.2  Turbulence  Measurement  on  a  Tornado 

There  is  a  current  flight  programme  on  the  Tornado  at  RAE  Bedford  for  the  measurement  of  structural  loads  and 
turbulence  response,  including  the  extraction  of  the  gust  time  histories.  The  structural  significance  of  combined  gust  and 
manoeuvre  loads  is  to  be  assessed  for  low  altitude  flight.  Comparison  of  the  flight  data  with  mathematical  model  predictions 
is  planned,  to  aid  the  understanding  of  the  combined  loading  actions. 


4.3  Helicopter  Measurements 

Under  this  section  on  low  level  turbulence,  reference  is  made  back  to  paragraph  2.5  where  a  helicopter  data  recording 
programme  is  discussed. 

5.0  CONCLUDING  REMARKS  AND  ACKNOWLEDGEMENT 
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